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Abstract 

We continue the study of the effective action for low x physics based on a Wilson renor¬ 
malization group approach. We express the full nonlinear renormalization group equation 
in terms of the average value and the average fluctuation of extra color charge density 
generated by integrating out gluons with intermediate values of x. This form clearly 
exhibits the nature of the phenomena driving the evolution and should serve as the basis 
of the analysis of saturation effects at high gluon density at small x. 
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Understanding of hadronic systems at large energies is one of the most exciting open 
problems in QCD today (see, e.g., a recent review J^). The recent wave of interest was 
sparked by experimental data on deep inelastic scattering [Q] showing a steep rise in the 
cross section at small values of Bjorken x. Theoretically it is believed that this rise is 
mainly triggered by fast growth of gluon density in a proton and indeed most hts to the 
HERA data use gluonic distributions that grow as a power at small x It argued a 
long time ago 0 that eventually the rate of growth of the gluon density should slow down 
and eventually saturate, thus curing a potential conflict with unitarity of the underlying 
scattering. Physically this should happen as the gluonic system becomes dense and 
gluons start overlapping in space causing self interactions, which are normally neglected 
in a linear QCD evolution, to become important. The most exciting problem of low x 
physics therefore is the question of how to deal with hnite density partonic systems. 

A formulation which is naturally directed towards this type of problems was suggested 
by McLerran and Venugopalan ^ in the context of large nuclei. Their approach was 
considerably modihed and generalized later §. [0 and evolved into an effective action 
approach to the low x hadronic structure. In this paper we continue to study this 
effective action with the help of Wilson renormalization group. The general Wilsonian 
renormalization group procedure for study of the low x physics was discussed hrst in p, 
and then in more detail in |^. It was also shown in the latter paper that in the linear 
regime, i.e. for small partonic densities, this renormalization group equation reduces to 
the celebrated BFKL equation for the evolution of unintegrated gluon density [0. This 
is a very important test for the consistency of the whole approach, however our goal is to 
be able to deal with the nonlinear regime where the color charge densities are not small 
and therefore have to be treated nonperturbatively. 

In this note we take one further step in this direction. The renormalization group 
procedure involves functional integration over two sets of variables: the fluctuation modes 
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of the gluon field with longitudinal momenta above the cutoff and the surface color charge 
density p at hxed value of p' which then becomes the charge density in the effective 
theory that describes only low longitudinal momentum glue modes. The goal of this 
note is to show how to perform one part of this program, namely the integration over 
the color charge density distribution. This will result in a gratifyingly compact and 
simple repesentation of the full nonlinear RG equation in terms of the average value and 
the average fluctuation of extra color charge density generated by integrating out gluons 
with intermediate values of x. The evolution towards smaller x is understood therefore 
in a simple way in terms of these two physical quantities. The simple form of the result 
is ideally suited to discuss the validity and the limitations of the BFKL equation from a 
new vantage point, as the low densitiy limit of the full nolinear RG equations. 

First, let us recall the framework as set up in [^. The starting point is the following 
action given in the Light Gone gauge = 0 

S = i J d‘^XtF[p'^{xt)] (1) 

- J d'^x^G"^ + ^ J d'^Xtdx~S{x~)p'"{xt)tTTaW_oc,oc[A~]{x~,Xt) 

Here is the gluon held strength tensor 

cr = di-Al - a‘-A^ + gUASA- (2) 


Ta are the SU{N) color matrices in the adjoint representation and W is the path ordered 
exponential along the x^ direction in the adjoint representation of the SU{Nc) group 


fF_oo,oo[^ ]{x ,Xt)=Pexp 


-ig 


dx'^A^ {x ,Xt)Ta 


( 3 ) 


The average of any gluonic operator 0{A) in the hadron is calculated as 

^ /T>[p°, 24^]0(24) exp{zg[p,24]} 

/T>[pGH^]exp{i^[p, A]} 


( 4 ) 
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The exponential of the imaginary part of the action 


lmS = J d^XtF[p'^] ( 5 ) 

can be thought of as a kind of “free energy” for the ensemble of the color charge density 
p{xt). The “Boltzmann factor” 

exp (^- J d^XtF[p^]^ (6) 

appearing in Eq.(^ controls the statistical weight of a particular conhguration of the 
two dimensional color charge density p°‘{xt) inside the hadron. 

In the classical approximation to the path integral in Eq.(^[^ one hnds the classical 
solution for the equations of motion at hxed p , and then averages over the charge density 
distribution with the “Boltzmann weight” Eq.(^. The classical solution for any hxed 
p{xt) has the structure]^ 

= 0 

= e{x-)a\xt) ( 7 ) 

where is related to p through 

Oii{xt) = {xt)diU\xt) 

ditti = -gp (8) 


and U{xt) is a unitary matrix. 

The quantum corrections to this classical approximation are large at small longitudi¬ 
nal momenta 0. To resum these large corrections we follow the renormalization group 
procedure as developed in 0 . 0 - 

^To be more precise, there is an infinite number of solutions to the classical equations of motion at 
fixed p. This is a consequence of the residual gauge invariance of the action. In the following we will 
be working in the gauge diAi{x~ —> —c») = 0. This is a complete gauge fixing and it therefore picks a 
unique solution of the equations of motion. 
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Let us introduce the following decomposition of the gauge field: 

A^x) = b;(x) + iA;{x) + a;(x) (9) 

where b'^{x) is the solution of the classical equations of motion, SA'^{x) is the fluctuation 
held containing longitudinal momentum modes such that < q~^ < Pn-i while 
is a soft held with momenta k~^ < P^, with respect to which the ehective action 
is computed. The initial path integral is formulated with the longitudinal momentum 
cutoh on the held 5A, < P^-i- The ehective action for is calculated by integrating 

over the huctuations 6A. This integration is performed with the assumption that the 
huctuations are small as compared to the classical helds More quantitatively, at each 

p+ p+ 

step of this RG procedure the scale is chosen such that In > 1, but as In 1. 

Expanding the action around the classical solution 6“(a;) and keeping terms of the hrst 
and second order in (5R we get 

^ = -^G{af - ^5A^[B-\p)r5A, + ga-p' + 0{{a-f) + iF[p] (10) 

where 

p' = p + 6pi + 6p2 (11) 

with 


Spl{xt,x^) 

and 


-2f^atSA‘i{x- = 0) 

Jdy* e(y* 


x'^) — 6{x'^ — y'^) SA '^{y~^,Xt,x 


( 12 ) 

0 ) 


Sp^,{x) = -r’’^[d^5A’l{x)]5A^{x) 

- j^p^ixt) J dy^SA~'^{y^,Xt,x~ = ^) J dz^SA~‘^{z"^,Xt,x~ = 0 ) 

X e{z^ - y+)e{y+ - a;+)trT“T"T'='T^ 

+ e{x+ - z+)e(z+ - |/+)trT“T^T"T^ 

+ e{z+ - x+)e{x+ - |/+)trT“T‘^T^T" (13) 
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The first term in both Sp1 and Sp 2 is coming from expansion of in the action while the 
rest of the terms proportional to p are from the expansion of the Wilson line term. The 
three terms correspond to different time ordering of the helds. Since the longitudinal 
momentum of a~ is much lower than of <5^4, we have only kept the eikonal coupling (the 
coupling to a~ component of the soft held), which gives the leading contribution in this 
kinematics. The inverse propagator is given by 


D-i 
D-i 
D-^ 


o 

ab 

- i+ 

. ab 
++ 
J ab 


x,y) 

x,y) 

x,y) 


+ D\b)D^{b)]^J^'^\x,y) 


(14) 


+ ‘^fabcal{xt)S{x )6{y )6^^\xt,yt)6{x^,y^)] 
{d^f5ab5^‘^\x,y) + fabcP%Xt)S{x~)5{y~)e{x'^ - y^)5^'^\xt,yt) 


The procedure now is the following: 


1. Integrate over SA'^ at hxed p and hxed 5p. 

2. Integrate over p at hxed p' = p + 5p. 


This generates the new ehective action which formally can be written as 

exp {iS[p\ a^^]) = exp - ^G^(a) + igap'^ (15) 

with 


exp {-F'[p']) = J D[p, 5A] 5{p' - p - 5p[5A]) exp 

Of course, to leading order in In l/x only terms linear in a* In l/x should be kept in F'. 
Dehning 


{-F\p]-'-5AD-^[p]SA) ( 16 ) 


O', 111 1a[p] = F'\p]- F[p] 

X 

gives the RG equation 

|;f|p]=a.A|p] 


(17) 


( 18 ) 
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In order for Eq. ([T8|) to be a bona fide renormalization gronp equation, we should be 
able to find the functional A for arbitrary Boltzmann weight F. That requires being 
able to integrate over the charge density p{xt) at arbitrary F[p]. The aim of this paper is 
to show how this is done. Let us start by considering the functional integral over SA at 
hxed p and p'. Although this integration is complicated |^, the structure of the result 
is simple and can be understood by simple counting of powers of the coupling constant. 
/^From the explicit expression of 5p in terms of SA, Eqs.(|T^), (0), it is readily seen that 


< Sp >SA= 0{as) 

< SpSp >SA= 0{as) (19) 


while all other (connected) correlation functions of Sp are higher order in Since we are 
working to the lowest order in we can neglect all these other terms. It therefore follows 
that to lowest order, after integrating over SA we are left with the weight function for Sp, 
which generates only connected one- and two-point functions. Such weight is obviously 
a Gaussian. Introducing the following notations 

1] := asln(l/a;) 

< Sp'^iXt) >SA =■ 

< 5p“(a;t,a;+)(5/(?/t,a;+) >5A =: vx''\xt,yt) (20) 


we can therefore write the result of the SAn integration in the form 


D[p, p'][Det{x)] ^/^exp(-F[p]) 


X exp [p(, -px- Wx] [Xxl] 

[ D[p, p'] exp{-t/[p,p']} 


Py Py V^y 


( 21 ) 


In the above equation we adopted condensed notations: the indices x stand for the set 
of indices and coordinates {xt,a}, and repeated indices are understood to be summed 
(integrated) over. We will use these notations in the rest of the paper. 
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We note here that this result can be derived formally by introducing the variable fl 


with the help of Lagrange multiplier 


(5(p' - p - (5p|(5/l]) = J £>1A] 


( 22 ) 


and subsequently integrating out A in perturbation theory to order Ug- 


The calculation of a and y is technically nontrivial |I3] and is beyond the scope of 
this note. Here we want to show, that assuming those quantities are known functionals 
of p, the integral over p in Eq.(^) can be performed explicitly. 

Consider again Eq.(^). We will assume that F[p] is a smooth functional of p and 
that the scale of its variation is independent of the strong coupling constant The 
last factor in Eq.(^T]) however is a very steep function of p which is sharply peaked 
around p' + 0{as). This is due to the factor l/a^ in the Gaussian weight for 6p. It 
is therefore clear that the integral can be calculated straightforwardly in the steepest 
descent approximation. We should also remember that we need to retain the terms up 
to hrst order in the coupling constant. Simple counting of the powers of a* shows that to 
reach this accuracy one needs to take into account up to the third order contributions in 
the steepest descent integration. In other words, the expansion of the exponential factor 
in Eq. (|^) is needed up to the fourth order in fluctuations around the stationary point. 


The steepest descent equation to hrst order in reads 

SF 1 _,6x 

^pu ^ ^pu 


Px px V^X VXst 


(23) 


Substituting this into Eq.( 
becomes 


we 


hnd that the logarithm of the integrand in Eq.(^ID 


U = F+-trln(y) 


(24) 


SF 1 , ,Sx' 

Opu ^ ^pu 


Xv 


SF 1 ^ ,6x, 

0 Py 2 0 Py 












In the above expression all the functionals are taken at - the solution of the steepest 
descent equation Eq. (p3D . 


Next we have to evaluate the correction due to Gaussian fluctuations of p around 
After some algebra we hnd 


5^U 

Py 


— ~'X ^ + 

A.XV ' 




p 

+x. 


^px^ P\ 


y L 


F+ -trln(x) 


(25) 


— I S(T', 


+ 


5py 5p, 


>U _l 

~Xuy 


+ 


— l^Xuv —l^Xuv 
Xxu -T— + X, 


6F 


-1 Sx 


5py 6p, 

Here the argument of all the functionals is again p^. 




Sp, 


S Pu 


As noted earlier to order we have to retain two more corrections to the steepest 


descent result. This is because the third as well as all higher derivatives of U in Eq.(^) 
is of the order l/a^ while the fluctuation is of order . Schematically, we can therefore 
write all order as contributions as 


D[p]exp{-U[p]) (26) 

= exp (-f/|p"]) / flip] exp [-\u"(p - py - - p»)’ - i(/""(p - p»)‘) 


= exp 


(-flip"]) /flip] exp (-Ifl"(p - py) 

1 - ^fl""(p - py +i (jy'"(p - p»)")' 


= exp 


(^-U[po] - ^trlnr/f/" - ^ (3 (jj) + (q (X^ + 6 O)) (27) 


Here we have used the obvious simplihed notations and graphical representation in terms 
of the “propagator” [U"]~^ and three - and four - point vertices U'" and U"". We have 
also re-exponentiated the result of the integration which is valid to leading order in 


Our task is simplihed, however by the fact that we need to retain only the leading 
contribution to this term. That for U'" and U"" we only need to keep the pieces of order 
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1/as, and for < (p — p°)” > only pieces of order Some straightforward algebra 

yields 


F. : = 


3 „ 3 6^U 

'6^U' 

-1 

' 5^U' 

-1 

4! 4! 5p„(5p„(5p^(5p2 

6p6p 

uv 

5p5p 

wz 


= h 


1. ( -1 X 1. i -i^x _i (5Xx 

T— ^)Xuv + -tr(x —X 


^ pu^ Pv 


^Pu ^pv 


1 ^ X.UV UW — 1 2 : 1 ; , 1<^X utt; — 1 Sx 


2 SpuSpy 2 ^p, 




6pv 2 (5p- 


Xwz 

V 


Spy 


(28) 


and 


K := 



2 3 ! ^pu^Pv^Pw ^Px^Py^Pz 


X 9 


6p6p 

+6 


-1 


(5p(5p 


-1 


fu_ 

6p6p 


-1 


yz 


'6^U' 

-1 

's^u' 

-1 

'6^U' 

1 

6p6p 

ux 

6p6p 

vy 

Sp6p 

WZ / 


= -V 


1. f -I 5x. w -i<^Xx ,1, / -i^X -i^Xx X 

otr(x ^JXuvt^ix ^J + Tt’^lX ^X ^)Xuv} 

^ S Pu S Py 4 0 Pii O Py 

1, , _1 5x X<5x uv l^x uv — 1 <^X , l^X uv — 1 ^x 

2 op^j OPj; 2 op^ ^Pz 2 OP;^ ^ Pu 


(29) 


Putting Eqs. (|25|) , (|26|) , (pq) and (p9|) together and eliminating the stationary value 
of p in favor of p' through Eg. (|23|) , we hnd 


F' = t/ + ltrln/^ + F3 + F4 

Px S Py 


= F 


Us ln(l/a;) 
2 

SF 6xuu 


r 

Xuv r e 
UpuSpv 




6F 6F 


C C C ')CuV r 

^ Pu S Pv S Py 0 Py 


wu. WA7/1) SCTy SF 

^2 ^ + 2^-2—(Ju 

SpU Spy Spy Spy 


(30) 


Equation (^) gives the Wilson renormalization group equation. 


d 


d\n.{l/x) 


F = — 


6F 6F 


p ^^Xuv 

Xuv C C r c c \uv r- 

SPuSPv SpySpy S Py S Py 




SpU Spy 

10 


spy 


^SF 

2— (Ty 

SpU 


( 31 ) 


























































This equation is extremely simple when written for the weight function Z = exp{—F} 

d 




(32) 


Equations (^ll) and are the central result of this paper. They provide the closed 
form of the renormalization group equation in terms of the functionals (j[p] and x[p]- 

We want to make now several remarks. First, the two terms in Eqs. (|3^) have a natural 
interpretation as the real and the virtual contributions to the evolution. The average 
fluctuation y is generated by the real ladder - like diagrams of the gluon fluctuation 6A 
while the average color density comes from the virtual one loop diagrams . 

Second, Eq.(^^ can be written directly as evolution equation for the correlators of 
the charge density. Multiplying Eq. (|3^) by Pxi---Pxn and integrating over p yields 
d 


dln{l/x] 


= a. 


< Pxi---Px„ > (33) 


^ < Px^--pXrn.,Pxrn+,-pX^_,Pxk + ,-pX„XxmX^ > 

0<m</c<n+l 


+ < Pxi---Pxi_iPxi + i---Pxn'^Xl > 

0<i<n+l 

In particular, taking n = 2 we obtain the evolution equation for the two point function 
d 


^ PxPy ^ Xxj/ T Px^y T Py^x 


(34) 


d\n.{l/x) 

In general this is not a closed equation since a and y depend on p in a complicated 
nonlinear way. However in the weak held limit, a and y are linear and quadratic func¬ 
tionals of p respectively. In this case Eq. (|3^ becomes a closed equation for the two 
point correlation function of p. The explicit expressions for a and y in this limit have 
been calculated in where it is also shown that Eq. (|3^ is precisely the famous BFKL 
equation. 

Since the correlator of the colour charge density is directly related to the unintegrated 
gluon density in a hadron (nucleus) P, Eq. can be straightforwardly rewritten as 
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a nonlinear evolution equation for the gluon density. It would be very interesting to 
compare this equation to the generalized evolution equation derived by Levin and Laenen 


rT|. In the second order in the gluon density (and neglecting the virtual contributions 


in the double logarithmic approximation) Eq. ( P^) should reduce to the famous GLR 
equation with the coefficient calculated by Mueller and Qiu . 


Another remark we want to make is that a general Gaussian form of the weight 
function 

F - j dx4vAxdi^-\x,,vM-m) (35) 

is not an eigenfunction of equation (^). This form of F was argued in ^ to be a good 
approximation for valence partons in a large enough nucleus, due to incoherence of color 
charges of partons coming from different nucleons. In general it should also be a good 
approximation in perturbative regime where F can be expanded in powers of p. We see 
however, that at low enough x non quadratic terms are generated through the evolution 
Eq. (pTD . In particular, even if we take the weak held expressions for a and y, the term 
quadratic in F on the r.h.s. of Eq.(PTD generates a correction quartic in p. Obviously, 
at weak coupling and small p this term can be neglected (or rather approximated by a 
quadratic term in a mean held like manner) so that a Gaussian Z is recovered in the 
perturbative regime. This suggests that a study of the RG how Eq.(pll) should indeed 
use Eq.(|3^) as an initial condition for the evolution starting at not too small values of 
X. The natural choice for the initial value of p{xt, yt) is 

r rt'^h 

pi{x„y,) = S{b,) j ^ m,Xo) (36) 

Here is the impact parameter, S(bt) is a nucleon shape function, xq is the value 

of X from which we start evolving F according to Eq.|^ and 4>{kt, x) is the unintegrated 
gluon density. 


Finally, we note that a and y are indeed calculable in a closed form. The result of 
this calculation will be presented elsewhere [ITU . 
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